function Fig2_g()
Fcolor='k';
Bcolor='w';
Fsize=18;  %font size

% %
RPA_k_25 = readmatrix('processed data/df_RPA_k_25.csv');
klb_25 = RPA_k_25(2:end,2);
phi_t_25 = RPA_k_25(2:end,3);
phi_25 = RPA_k_25(2:end,4);

RPA_k_70 = readmatrix('processed data/df_RPA_k_70.csv');
klb_70 = RPA_k_70(2:end,2);
phi_t_70 = RPA_k_70(2:end,3);
phi_70 = RPA_k_70(2:end,4);

f = figure();
f.Color = Bcolor;
hold on
plot(klb_25, phi_25)
plot(klb_25, phi_t_25)
plot(klb_70, phi_70)
plot(klb_70, phi_t_70)
ax = gca;
box on
ax.Color=Bcolor;
ax.XColor=Fcolor;
ax.YColor=Fcolor;
xl=xlabel('x');
yl=ylabel('y');
ax.FontSize=Fsize;
yl.FontSize=Fsize;
xl.FontSize=Fsize;
axis([0.25,5,-2.5,0.1])

%%
ED_r_25 = readmatrix('processed data/df_ED_r_25.csv');
rs = ED_r_25(2:end,2);
phi_2deg = ED_r_25(2:end,3);
phi_imp = ED_r_25(2:end,4);
f = figure();
f.Color = Bcolor;
hold on
plot(rs, -phi_imp+phi_2deg,'LineWidth',0.5,'DisplayName','Vtotal')
plot(rs, -phi_imp,'LineWidth',0.5,'DisplayName','Vimp')
plot(rs, -phi_2deg,'LineWidth',0.5,'DisplayName','-V2deg')
legend();
title('V2deg')
ax = gca;
box on
ax.Color=Bcolor;
ax.XColor=Fcolor;
ax.YColor=Fcolor;
xl=xlabel('r (nm)');
yl=ylabel('phi (meV)');
ax.FontSize=Fsize;
yl.FontSize=Fsize;
xl.FontSize=Fsize;
axis([0,150,-20,5])



end